Small amplitude quasi-breathers and oscillons 



O 
00 

o 
X 



Gyula Fodor^, Peter Forgacs^'^, Zalan Horvath^, Arpad Lukacs^ 
^MTA RMKI, H-1525 Budapest 114, P.O.Box 49, Hungary, 
^Institute for Theoretical Physics, Edtvos University, 
H-1117 Budapest, Pdzmdny Peter setdny 1/A, Hungary, 
^LMPT, CNRS-UMR 6083, Umversite de Tours, Pare de Grandmont, 37200 Tours, France 

Quasi-breathers (QB) are time-periodic solutions with weak spatial localization introduced in 
G. Fodor et al. in Phys. Rev. D. 74, 124003 (2006). QB's provide a simple description of oscillons 
(very long-living spatially localized time dependent solutions). The small amplitude limit of QB's 
is worked out in a large class of scalar theories with a general self-interaction potential, in D spatial 
dimensions. It is shown that the problem of small amplitude QB's is reduced to a universal elliptic 
partial differential equation. It is also found that there is the critical dimension, Dcrit = 4, above 
which no small amplitude QB's exist. The QB's obtained this way are shown to provide very good 
initial data for oscillons. Thus these QB's provide the solution of the complicated, nonlinear time 
dependent problem of small amplitude oscillons in scalar theories. 
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I. INTRODUCTION 



Static scalar lumps (finite energy particle-like solutions in scalar field theories) are known to be absent in more than 
two spatial dimensions (D > 2), however very long lived oscillating lumps have been observed in scalar theories with 
rather general self- interaction potentials in spatial dimensions D < 7 [H, 0, d, 0, Si ■ These states (baptized oscillons in 
Q^' Ref. 3) are of quite some interest, in spite of the fact that they eventually decay, since oscillons evolve from rather 
D [ generic initial data in a remarkably large class of theories. 

In 2 spatial dimensions, extremely long living breather-type objects have been found in the sine-Gordon (sG) model 
I 0, and more recently in Ref. Q generic oscillons, with lifetimes of the order of 10^ in natural time units have been 
^— I . exhibited both in the sG and in a (j)^ scalar model. In 3 spatial dimensions spherically symmetric oscillons in the (j)'^ 
J> ■ scalar model with lifetimes « 10'' have been found and investigated in detail [1, [^, [lO'|. Recently oscillons have been 
numerically observed starting from rather generic initial data in spontaneously broken SU(2) gauge theories (for the 
■ case when the mass ratio of the Higgs field to the W boson is « 2), [11]. More recently 3 dimensional numerical 
[ simulations have shown [T5| that long living oscillons are present in the SU(2)xU(l) model, i.e. in the full bosonic 
. ' sector of the Standard model of electroweak interactions, which is clearly of great potential interest. According to the 
findings of Ref. [l3| non-symmetric oscillons evolve towards symmetric ones (at least in 1 -I- 2 dimensions) , indicating 
that the long time evolutions will be dominated by spherically symmetric configurations. 

Oscillons are already interesting in their own right as very long living lumps, appearing in many physical theories, 
whose existence is due to the non-linearities. They are expected to have important effects on the dynamics of 
various systems (including the Early Universe), since they retain a considerable amount of energy. Importantly, 
oscillons can easily form in physical processes such as the QCD phase transition, where oscillon like objects in 
the axion field have been observed [l4|, in vortex-antivortex annihilation and in domain collapse [l^. In [l7| 
^ oscillating field configurations (I-balls) were found in potentials where the quadratic term dominates. Once formed, 
oscillons could considerably influence the dynamics of the system as has been suggested for the case of the bubble 
nucleation process p^ . Oscillon formation has been reported after supersymmetric hybrid inflation fl^. A slightly 
different mechanism for the formation of long lived objects (quasi-lumps) during first order phase transitions has been 
investigated in |20| . The persistence of oscillons in one spatial dimension in an expanding background metric has been 
reported in Ref. |2l|. M ost recently in a study of semiclassical decay of topological defects possible oscillons formation 
has been reported |22| . 

Infinitely long lived oscillons with finite energy, commonly known under the name of breathers, are rather excep- 
tional. Simple heuristic arguments indicate that spatially localized time periodic solutions (breathers) do not exist 
in generic theories [1^, and that under a general perturbation breathers are unstable [2^. In 1 spatial dimension 
the absence of small amplitude breathers in a scalar theory with (f>^ interaction has been demonstrated in Refs. [25| 
[2^ . More generally it has been shown in Ref. ^] that the only non-trivial 1-1-1 dimensional scalar theory with real 
analytic self-interaction potential which admits breather solutions is the sine-Gordon (sG) theory. Interestingly in 
a recent work [28^ a whole family non-radiating solutions of the 1 + 1 dimensional 'signum-Gordon' model has been 
found. This illustrates that one can expect new and surprising phenomena in theories where the interaction potential 
is not smooth. 
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On the other hand, bounded, time periodic solutions of nonUnear wave equations in R are abundant. By per- 
forming a Fourier decomposition in time, the problem of finding bounded, time periodic solutions is reduced to solve 
an infinite set of coupled nonlinear partial differential equations (PDE's). Even the simplest cases, such as = 1 or 
spherical symmetry, when the Fourier amplitudes satisfy an infinite set of ordinary differential equations (ODE's), are 
non-trivial to analyze, but it is clear that a plethora of bounded solutions exist. The class of bounded, time periodic 
solutions in R-^ contains families of breather-like objects, which have a well defined core, outside of which the fields 
fall of rapidly, but barring exceptional cases, in the far field region there is also a radiative tail, which corresponds to 
a standing wave. Because of the asymptotically standing wave asymptotics, such objects are only weakly localized in 
space. The energy density of the core is much larger than that of the tail, however, the total energy contained in the 
tail is infinite. Considering such a weakly localized object in a finite volume, V, its total energy, is proportional 
to E (X V^l^ . Such objects can be thought of as radiating lumps, whose energy loss is compensated by a flux of 
radiation coming from infinity, rendering the system time periodic . In [lo| a special class has been singled out in 
the huge phase space of bounded, time periodic solutions, obtained by minimizing the energy in the radiative tail. 
Solutions belonging to this class have been called as quasi-breathers (QB's). Even if the QB's are not physical objects 
by themselves in the whole space, R^, because of containing an infinite amount of energy, nevertheless a finite piece 
of them (containing the core and even part of the tail) constitutes a very good approximation to oscillons in a R^, 
as demonstrated in detail in the scalar theory in D = 3 10]. 

In this paper we carry out a systematic analysis of bounded, time dependent solutions in a -D dimensional scalar 
theory with a rather general class of self-interaction potential in the limit when the amplitude, e, of the solution goes to 
zero. This way we obtain a rather general method to find small amplitude QB's. Following Refs. [2^,[2^[23] we derive a 
formal series solution in e whose terms are all bounded both in time and space. We show that demanding boundedness 
in time necessarily leads to periodicity, and that the time and space dependence of the solution separates. We derive 
a single master equation determining the spatial dependence of the leading term in the series, which turns out to be 
universal for the class of scalar theories we consider. This equation is a nonlinear elliptic PDE with a cubic nonlinearity. 
It turns out that exponentially localized solutions of the master equation exist in spatial dimension, 13 < 4. While in 
D = 1 the solution is unique, for D > \ there is an infinite family of exponentially localized solutions. In the case 
of spherical symmetry, members of this family can be characterized by the number of their nodes. Solutions with 
nodes contain considerably higher energies than the fundamental one, nevertheless they also correspond to oscillons. 
The higher order terms in the small amplitude series can be obtained from linear inhomogeneous PDE's whose source 
terms are determined by the localized solution of the master equation. Each term of the small amplitude expansion 
obtained this way is exponentially localized, and one could think that it represents a breather. In general this series 
is, however, not convergent, it is rather an asymptotic one. This fact refiects the absence of 'genuine' breathers 
with spatial localization. One can think of this series as an excellent approximation of QB's whose radiative tail is 
smaller than any power of £. Therefore the small amplitude series provides only the exponentially localized core part 
of the QB. Since the amplitude of the radiative tail is so small this does not really matter and in fact the oscillon 
states corresponding to such initial data for small values of e do have very long lifetimes. Even the first few terms 
in the small amplitude expansion yield quite good initial data for long living oscillon states, as demonstrated by our 
numerical simulations for the standard 0* theory for the spherically symmetric case in Z3 = 2 and 13 = 3. As already 
mentioned, one of the interests of the small amplitude QB's, is that they can be identified with the core part of small 
amplitude oscillons, which radiate very weakly, and hence they have a very long lifetime. We have verified this by 
numerical simulations, namely we shown in dimensions D = 1 , 2 , 3 that small amplitude QB's do provide excellent 
initial data for long lived oscillons. 

We have also computed the energy of the QB's approximated by the small amplitude series, which corresponds to 
the energy content of their core and this is of course always finite. We have shown that the energy of the QB core 
is a monotonously decreasing function of the frequency near the mass threshold in dimensions D < 2. This implies 
the absence of a critical frequency minimizing the energy in dimensions _D < 2, in agreement with known numerical 
results. For D > 2 the energy of the small amplitude QB's increases without bound as their frequency approaches 
the mass threshold. This implies the existence of a critical frequency where the energy is minimized. 

As it has been already mentioned, well (exponentially) localized QB's with a finite energy core exist only for D < 4. 
In higher dimensions, D > 4 the exponential localization property of the core of the QB's is lost and in fact the very 
existence of a well defined core is problematic. In Ref. fHl it has been already suggested that oscillons cease to exist, 
for dimensions greater than 5 or 6 depending on the details of the potential. According to the findings of Ref. [2^ 
the lifetimes of oscillons decreases rapidly as the dimensionality of space is increased, and the QB picture does not 
give a good description. We have pushed further the analysis of small amplitude oscillons in higher dimensions to 
understand the situation better. We shall present our results in dimensions D > 4 in a sequel to this paper [30| . 
Without going into too much details we can state the following. Small amplitude oscillons do exist in dimensions 
Z3 > 4, without any apparent limitation for D. These small amplitude oscillons in dimensions Z) > 4 are, however, 
qualitatively different from their lower dimensional counterparts. In particular, they are not well (exponentially) 



3 



localized, and they cannot be described by QB's in the sense of Ref. flO[. The energy of higher dimensional smaU 
amphtude osciUons also becomes quite large, therefore they are probably less interesting physically than genuine QB's, 
however, by choosing suitable initial data, they can still have very long lifetimes. Their large energy content may 
explain that such osciUons in > 4 are somewhat more difhcult to be found. 

The plan of our paper is the following: In Section [IT] first the class of scalar models to be studied is introduced, 
then the small amplitude expansion in D spatial dimension is carried out. We derive the master equation without 
any symmetry assumptions, and calculate the QB solution to order 4 in the e expansion as well. The energy of the 
QB's is also computed, and in the last subsection the existence of the critical dimension D = 4 is derived above 
which no small amplitude QB's exist. Section Hill is devoted to a detailed numerical analysis of the solutions of the 
spherically symmetric master equation and the explicit computations of the higher order terms in the e expansion in 
spatial dimensions D = 2 ,3. In Section HVl results on the numerical time evolution of QB initial data up 6th order in 
e is presented for D = 2 ,3. 



II. THE SMALL AMPLITUDE EXPANSION 



In this Section we carry out a detailed analysis of the small amplitude limit of QB's of the NLWE ^ in D spatial 
dimensions, without any symmetry assumptions. It turns out that bounded non-trivial small amplitude solutions 
are periodic in time, and that the time and spatial dependence completely separates. In subsection III Bl we derive a 
universal elliptic PDE, referred to as the master eq. governing the behaviour of the solutions. Next, in subsection 
III CI the solution is obtained up to order 4 in e in a general class potentials, and up to order 6 for theories with a 
symmetric potential, such as the sine-Gordon model. 



A. The class of theories considered 

We consider a scalar theory in a 1 4- _D dimensional flat Minkowski space-time, with a general self-interaction 
potential, whose action can be written as 



A = dtd 



(1) 



where (/) is a real scalar field, dt = d/dt, di — d/dx^ and i = 1, 2, . . . , I?. The equation of motion following from ^ is 
a non-linear wave equation (NLWE) which is given as 

-^,tt + A^^U'i<f>)=^ + Y,9k^''. where A = ^ . (2) 

k=2 i=l "^^ 

In Eq. (21) the mass of the field is chosen to be 1, and it has been assumed that the potential, [/(</>), can be written as 
a power series in 0, where the gu are real constants. For the standard 0^ theory this interaction potential is simply 

U{<t>)^\(p^{<t>-2f , U'{(p)^(p-^-4>^ + \^\ (3) 

i.e. 52 = — |, 53 = ^ and 5^ = for i > 4 in this case. Note that in our previous paper, Ref. a different scaling of 
the (f)^ potential has been used, making the value of the mass to be m = y/2 instead of the present value m = 1 used 
in this paper. For the sine-Gordon potential U{(j)) = 1 — cos(0), we have g^i = and 521+1 = (~1)V(2« + !)!■ 
The energy corresponding to the action ([1]) can be written as 

E = j d^xE, f = 1 V i {d,4>f + C/(0) , (4) 

where 8- denotes the energy density. 

B. Derivation of the master equation 

We are looking for small amplitude solutions, therefore we expand the scalar field, 0, in terms of a parameter e as 

(^ = ^£Vn. (5) 

fe=l 
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In order to obtain non-trivial solutions of Eq. ([2]) their characteristic scale must also become e-dependent. The size 
of smooth configurations is expected to increase for decreasing values of e, therefore it is natural to introduce new 
spatial coordinates by the following rescaling 

C = ex' . (6) 

One must also allow for the e dependence of the time-scale of the configurations, therefore a new time coordinate is 
introduced as 

r = u;{e)t. (7) 
uj{e) is assumed to be analytic near the threshold, uj — I, and it is expanded as 



u;\e)^l + ^e''LUk. (8) 

k=l 

After these rescalings Eq. ([2]) takes the following form 

oo 

-uj^ + e^Acj^^cj^ + Y^gkcf^" ■ (9) 

In equation ^ and in the rest of this section an overdot stands for the derivative with respect to the rescaled time 
coordinate, r, and all spatial derivatives are taken with respect to the rescaled coordinates C. Substituting the e 
expansion of the scalar field, (f>, and of uj'^ into ^ the equations determining the first three lowest order terms are: 

4>i + (^i - 0, (10) 

(f'2 + (f>2 + 92Ct>l + ^l4>l = 0, (11) 
03 + 03 + 2g201</'2 +530? - 01 - A(?!)i + Wi02 + t^201 = 0. (12) 

As it is clear from Eqs. PTH) - P^ the time dependence has been separated from the spatial one, and we have obtained 
a set of harmonic oscillator equations. Now the solution of Eq. pH)) is clearly given by 

01 Pi cos(t + a) , (13) 

where pi and a are functions of the spatial variables C*- That is, the lowest order term of the solution is just a 
harmonic oscillator in time, with frequency uj = 1 (note that this is the same with respect to both time coordinates, 
r and t at this order). This distinguished value of the frequency, uj — 1, corresponds to the threshold determined by 
the mass of the scalar field. 

As Eq. (jlip is a linear inhomogeneous equation, its solution is easily obtained: 

02 = P2 cos(t + a) + q2 sin(r + a) + ^p\ [cos(2t + 2a) - 3] -|- [2t sin(T + a) + cos(t + a)] . (14) 

6 4 

Since we are looking for bounded solutions, it is necessary to impose uji — Q. This amounts to demanding the absence 
of the resonance term tJi0i in eq. (jlip . Substituting the solutions for 0i and 02 into Eq. (I12p one obtains yet another 
forced oscillator equation for the time dependence of 03: 



^3 



(piAa + 2VaVpi) sin(T + a) — [Api + UJ2P1 + Ap? — pi(Va)^] cos(t + a) 



+ ^P?(25l + 3<?3) cos(3r 3a) + sin(2T + 2a) + P2 cos(2t + 2a) P2] = , (15) 

where we have introduced the combination 

5 3 

A = gff2 - ^53 , (16) 

which will play an important role in the following. For the standard (f>^ theory this parameter takes the value 
A = 3/2. For the sine-Gordon potential A = 1/8. 

As already explained we are looking for bounded solutions, therefore it is necessary to guarantee the absence of 
resonance terms also in Eq. (fT5|) . The vanishing of the coefficient of the sin(T -I- a) terms implies 



V(p2Va) = 0, (17) 
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and from this equation one immediately derives the following condition 

/ aV(piVa) = / api^-Va- / Pi(Va)^=0. (18) 
Jn J on Jn 

Assuming that the integrand of the boundary term vanishes sufficiently fast we conclude that Va = 0. Our assumption 
is quite reasonable since we are looking for bounded solutions in time and localized in space. Therefore a must be a 
constant which can be absorbed by a shift in the time variable. From now on we set a = 0. Then the vanishing of 
the coefficient of the resonance term proportional to cos r implies 

Api + LU2P1 + Xpl^O. (19) 

A necessary condition that this equation admit exponentially localized solutions is a;2 < 0, which we shall assume 
from now on. In this case we can set a;2 = — 1 by a simultaneous rescaling of and pi. This rescaling corresponds 
to choosing a different parametrization e for a solution with a specific frequency uj. 

By an analytic redefinition of the expansion parameter, e, and by rescalings, all coefhcients LOi can be made to 
vanish for i > 2. This means that by a suitable transformation of the expansion parameter, e, one can always achieve 
that the following relation between the frequency, w, and the expansion parameter, e, holds: 

Lo^V^-e^- (20) 

In terms of the physical time coordinate, t, the configuration oscillates with frequency uj. This is the physically 
important frequency characterizing these periodic solutions. Apart from the leading order behaviour, the precise 
choice of how s depends on uj is physically irrelevant. 

After setting i02 = —1, it is easy to see (multiplying p9p with pi and integrating) that A > is a necessary condition 
for the existence of bounded solutions vanishing at infinity. Assuming A > 0, by rescaling pi one obtains 

AS* -5 + 5*^ = 0, 5 = pi\/A. (21) 

This master equation constitutes an important result equation of our paper. Quite remarkably Ea. (|2ip is universal 
for the class of theories considered, and the dependence on the parameters of the interaction potential enters only 
through the combination A when reconstructing . 



C. Higher orders in the e expansion 

Let us now turn to the determination of some higher order terms in the e expansion. Armed with the simple form 
of the solution for pi it is now easy to obtain the explicit time dependence of ^3 by integrating Eq. ([1| 



03 = gasinr + P3COST+ y |i (^gl - A^ cos(3t) +32 [92 sin(2T) +p2(cos(2r) - 3)]| . (22) 

In 03 two new functions, p^, qs have appeared. The absence of resonances at fourth order yields two conditions: 

Ag2P? - 92 + A(72 - , (23) 

which is the vanishing condition of coefficient of the term of sinr, and another one 

3Ap2P? - P2 + AP2 = , (24) 

which ensures the vanishing of the coefhcient of cost. 

Let us first note, that quite remarkably q2 oc pi actually solves Eq. (j23|) . Therefore by shifting the time coordinate 
by a suitable term of order e we can eliminate (72 ■ In higher orders of the e expansion, we have verified that up to 
order the vanishing of the coefRcient of the terms proportional to sinr leads to equations which are equivalent 
to Eq. ((23|) . We conjecture that this is in fact true to all orders, i.e. by a suitable choice of the origin of the time 
coordinate r one eliminate all terms proportional to sinr. This observation is quite important because it implies that 
all small amplitude QB-type solutions necessarily possess time reflection symmetry. This is of course what one would 
expect based on simple physical intuition. Let us point out here, that all long-lived oscillon configurations observed 
in time evolution simulations appear to show this symmetry to a very high degree. Of course oscillons are not exactly 
time reflection symmetric because they radiate some energy to infinity. Moreover, time reflection symmetry has been 
usually implicitly assumed when performing Fourier decomposition in order to flnd time periodic states. For all these 
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reasons it would be of interest to find a mathematical proof of the validity of time reflection symmetry for periodic 
solutions of NLWE's. 

There is no reason to expect that Eq. (|24p admits bounded solutions vanishing at infinity apart from those corre- 
sponding to the translational and rotational symmetries, therefore from now on we set p2 = 0. The vanishing of the 
coefficient of the cos t term in the fifth order equation yields 



9l 



Ap3 -P3 + 5XpIp3 + yPl {I7pl + 19(Vpi)2 



Pi 

216 



(3785452 + 36A5I - 28O5I - - 13555) = . 



(25) 



This is a linear, inhomogeneous equation. It can be brought to a much simpler form by introducing a new variable Z 
instead of p^: 



1 



P3 



24 



Then Eq. ((25)) takes the compact form 



-M 



5 7 

g55 - ^5254 



— a'^]Z 
27 ' 



193'^) 



AZ - Z + 3S^Z - S'^ 



0. 



For the specific example of the standard ip^ theory this relation is simply 



P3 



3V3 



(26) 



(27) 



(28) 



Integrating the corresponding equation for 05 a new unknown function, p^ appears, in analogy to the third order case. 

To summarize, we have obtained the solution of the NLWE ^ in the small amplitude expansion up to order four. 
All terms have harmonic time dependence, and the spatial part is determined by the two universal elliptic PDE's, 
eqs. (|21l27p . There is no obstacle to continue the computation to higher orders, the general formulae become then 
quite complicated of course. The small amplitude expansion of the solution of Eq. ([9]) up to order four for general 
interaction potentials can be written as: 



Pl COST 

^52P? (cos(2t)-3) 



1 

72 



(452-3A)p?cos(3r) 



6 

P3 cos T 

(354 - 552A + 55I) cos(4r) 
(852(Vpi)2 



Ugipl + Wglpj - 2452P1P3 



23g2Xpt - 852P?) cos(2r) 
3 



(29) 
(30) 

(31) 
(32) 



1 31 o 

~92pj - 52P1P3 + g52Apt - 52(Vpi)2 + -j^glpf - -^9iPt 

A considerable simplification occurs when the scalar self-interaction potential, U{(j}), is symmetric around its min- 
imum (f) = 0. In this case 524 = for all i = 1 , . . ., and all even power terms in the s expansion vanish, i.e. 02j = 
for i — 1 , . . .. Since 02n contains only terms of the form cos(2A:t) with k = 1 , . . .n, and (j)2n+i contains only terms 
proportional to cos((2fc -I- 1)t) , with k = 1 , . . .n, this also implies that for such symmetric potentials no even terms 
in the Fourier expansion arise. In this case ps is proportional to Z and the equation determining the function p^ 
becomes reasonably simple, it can be written as 



Ap5 - P5 



iS^P5 



SZ 
576VA 



(3Z - 



V A2 
/ 31557 6O55 



576VA V 



A2 



(33) 



Then 



= P5 COS T 

s 



55 



384\/A 



1152VA V -^^ 
3O55 
A2 



f% + 2)cos(5r) 



2 5Z+12S'^-12(VS')"- 



1555 
A2 



-2 5* 



cos(3r) . 



(34) 
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These expressions encompass for example the case of the sine-Gordon model. The corresponding equations for the cj)'^ 
theory, in the case of spherical symmetry, will be listed in Section Hill 

As we have already stressed several times, it is by now well understood that spatially localized breathers of the 
NLWE ^ do not exist in for general analytic potentials, even if a general mathematical proof is known only in 
D — 1. Let us note here, that a remarkable example in D = 1 admitting non-radiating breather-type solutions in the 
framework of "^''-shaped (non-differentiable) potentials evades this theorem [2^. In the case of analytic potentials, 
where the theorem applies there is still a point to be stressed. Assuming that exponentially decreasing solutions of 
the master equation exist, all higher order terms in the small amplitude expansion are also exponentially localized, 
and they are periodic in time. As we have learned from the example of the one dimensional (j)^ theory [2^ the series 
solution in powers of e does not converge to a breather, it is an asymptotic series. Nevertheless, to a given order 
in the expansion for sufficiently small values of e the corresponding sum yields a configuration with a spatially well 
localized core. This time periodic configuration corresponds to a QB whose standing wave tail is smaller than e" for 
any n > 0. As it will be shown in Section HVl such QB's constitute an excellent approximation to an oscillon. 

D. The energy 

In this subsection we evaluate the energy of small amplitude QB's, in D dimension. In the rcscalcd coordinate 
system, r , ^, the energy of a configuration, Eq. ([4]) can be written as 

E^^Jd^'CS, where £ ^ l{l ^ e^) {drcl>f + e^^ {d,cl)f + U{(l>) . (35) 

Because of the periodic time dependence we shall compute the energy density averaged over a period, 

1 r^^ 

E=— drE, (36) 

and in this subsection the bar over a quantity will denote its time average. Using the results of the e expansion, Eqs. 
([29|) - ((3T|). the time averaged energy density, up to fourth order in e can be written as 



£ = --S 



2 



XS'^i'o^gl + 27A)(S'2 2) - b^X^iySf 

(37) 



2A 216A3 

-5*^(13535 - 378.g2.94 + 280g| - m\gi + 9A^) 



For the (f)^ theory the time averaged energy density takes a much simpler form: 

. (38) 



\(ySf - —S\S^ + 2) + —SZ 
&^ ' 108 ^ ^36 



Using the above results, the time averaged total energy, E, has the following e expansion in D dimension: 

E = ^ ^i-Dj^^ ^ ^j^g^g ^ y g2 ^ (39) 

and El denotes the integral of the 4-th order term in the energy density, ([37]) in . As one sees from Eq. ([39|) 
the leading order behaviour of the time averaged total energy is E (x e^~^ . This implies that the e-dependence of 
E changes essentially at D = 2. In dimensions D > 2 the total energy increases without any bound for decreasing 
values of £. In D = 2 E tends to a constant, and for _D < 2 it goes to zero as e ^ 0. This also implies that the core 
energy of a QB in dimensions D > 2 should exhibit a minimum for some frequency uj^. In fact, from Eq. (j39p one 
immediately finds 

''^-^''^--^'2XiA~D)Ei- 

The above result can only be taken as an indication of the minimum even if Em 1. The numerical values of Eq and 
El will be given for the fundamental solutions in D = 2 and D = 3 in case of spherical symmetry in Section IIIII We 
note that in D = 1 S{C) = \/2sech(C), therefore Eq = 4, and the leading order e dependence of the energy is given as 
E = 2e/X + 0{e^). 
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E. Critical dimension D — 4 



In the following we present sonic simple, although important results concerning the existence of spatially localized 
solutions of the master equation, (HH). First spatially localized solutions of (PTjl which have a limit for ja?"! — ^ oo 
should decrease exponentially, since those tending to a constant exhibit oscillatory behaviour. Next we show that 
exponentially localized solutions of Eq. ((2T|) cannot exist for D > A, implying that small amplitude QB's exist only 
in dimensions D < 4. To see this, consider the following virial identity derived from equation (|2ip : 

((V5)2) + (52)-(54)=0, (41) 

where (/) :— J dP xf{x). Furthermore, another virial identity can be found from the scaling transformation {x fix) 
of the action corresponding to (EI]), / d^x[{VSf + S'^ - S^/2]: 

{D-2){(ySf) +D{S^) ~^{S^) ^0. (42) 

From Eqs. (|^T|) and ([1^ one immediately finds 

2{S^) + ^{D-4){S')=0, (43) 

which equality can only be satisfied if _D < 4. 

The absence of small amplitude QB's in more than 3 spatial dimensions does not imply per se that oscillons would 
be also absent if 13 > 4. As a matter of fact we have found that small amplitude oscillons exist in dimensions D > 4, 
without any apparent limitation on D. These higher dimensional {D > 4) small amplitude oscillons do not have a 
well defined, exponentially localized core, and they cannot be described by small amplitude QB's. Even the total 
energy of the core in > 4 is not well defined. Interestingly by choosing suitable initial data, for increasing energy 
content one can achieve that they have very long lifetimes. Various arguments and numerical studies of spherically 
symmetric oscillons in Z3-dimensions by Gleiser 5] led him to conjecture the existence of a critical value of D (Dc ^ 6) 
above which no long lived oscillon states would exist. The existence of higher dimensional small amplitude oscillons 
contradict this conjecture, however, since these contain very large amount of energies it may have been less obvious 
to start with such initial data. This might explain why such objects have been missed. Also for a fixed amount of 
energy, the lifetime of oscillons exhibits a significant decrease for D > 3. The results of the recent work [29] show 
that in _D = 5 the lifetimes becomes as small as a few 100 (in natural units). Let us mention here another interesting 
point. In Rcf. a very long lived oscillon state has been exhibited in D = 6. This object is not in the class of small 
amplitude oscillons, and should be understood better. In any case this provides another example how rich the phase 
space of time dependent solutions of a simple non-linear wave equation can be. 



III. SOLUTION OF THE MASTER EQUATION FOR SPHERICAL SYMMETRY 

In this Section we consider spherically symmetric configurations, in which case the PDE's determining the functions 
S, Z, etc. reduce to ODE's. This simplifies of course significantly the problem of solving both the master equation 
and the associated inhomogeneous ones. We exhibit some numerical solutions of these equations in 13 = 2 and in 
D — 3. We present the solution of the e expansion in the (p^ theory up to 6th order. 

For spherically symmetric configurations the master equation (|2ip takes the form 

d^S D-ldS ^ ^3 ^ 

TT + — -S + S^^O, (44) 

dp^ p dp 

where S" is a function of the rescaled radial coordinate p = er. In 1 spatial dimension the solution of (|44p vanishing 
at infinity is unique, it is given explicitly by = v^sechp. In contradistinction to D = 1, for higher dimensions, 
1 < D < 4, the solution vanishing at p — > oo is not unique. Our numerical analysis indicates, that for 1 < D < 4 
there is a family of localized solutions of Eq. (|44p indexed by the number of zeros (nodes) of S{p). On Figures [T] and 
[5] the first few members of this solution family are exhibited in dimensions D ~ 2 and D = 3. As it has been shown 
in the previous Section, there are no solutions of Eq. (|44p for D > 4, which tend to zero for p oo. The values of 
S at the origin p = are tabulated in Table IH It can be already expected that the fundamental solution (without 
nodes), Sq, is physically the most important. Indeed, as we shall show later, oscillons corresponding to solutions of 
the master equation with nodes contain more energy and have significantly smaller lifetimes than those corresponding 
to nodeless ones, and they are also less stable. 
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FIG. 1: Solutions of the master equation (|44p in D = 2 with 0, 1, 2 and 3 nodes 
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FIG. 2: Solutions of the master equation (|44|l in D = 3 dimensions with 0, 1, 2 nodes. 



number 
of nodes 


Sip = 0) 
D = 2 D = 3 





2.20620086 


4.33738768 


1 


3.33198927 


14.10358440 


2 


4.15009404 


29.13121158 


3 


4.82960282 


49.36070988 



TABLE I; Central values of S for two and three dimensions. 
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In the case of spherical symmetry, equation (|27p determining the third and fourth order terms in e takes the form 

(fZ D -1 dZ 



, , , Z + iS^Z - ^Q. (45) 

dp'^ p dp 

Since Eq. (j45p is Unear for the unknown, Z, with inhomogeneity it admits a globally regular solution for any S 
regular at p = and vanishing for p ^ oo. Some numerical solutions of Eq. (|i5|) are plotted on Figs. [3] and S] in two 
and three dimensions. 




012345678 



P 

FIG. 3: Solutions of Eq. (|45|) for Z corresponding to S without and with 1 node in two dimensions. The central values are 
1.45076 and 0.6018575, respectively. 



A. Solution for the (j)^ theory up to 6th order 

As already discussed, the equation determining the function, Z, is universal for any choice of the potential U{(j)). 
Of course the terms in the small amplitude expansion do depend on the potential U {(p) , and the reconstruction of 
^3 from Z using Eqs. ([26]) and ((3T|) depends of the values of the coefhcients gi. In the rest of this paper we shall 
concentrate on the 4>'^ theory given by Eq. ([3]) and provide the results of numerical simulations only for this case. 

The equation determining the fifth order term in the small amplitude expansion, cf)^^ through can be written as 

dp'^ p dp 64 

^53^^-1)^ fdSV ^ mn-l),.^,. _ ,^dS ^ 8287^, ^ ^ ^ ^^^^ 
\dpj p dp 

where Y is defined by 



\/2 f^^ 1235^2^ , 1503 ^ 17^3 11525 . 



S'Z+ ^-245* 5^ + . (47) 

9^/3 V 32 16 3 384 / ^ ^ 

The solution of Eq. corresponding to the fundamental solution of the master Eq. 5*0, has the central value 
Yo{0) = —87.78183 in two dimensions and ^0(0) = 60356.38 in D = 3. These actual values themselves have no 
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FIG. 4: Solutions of Eq. (gSJ in three dimensions. The central values of Z are -16.17403 and -1290.021. 



physical significance in view of the scaling freedom in the definition of the function Y . Clearly, any constant could 
have been included in front of the term Y in Eq. (|47|1 without modifying the final results for the magnitude of the 
(^i's in the e expansion. 

In the following we list the values of the terms of the small amplitude expansion, <^i, up to order six in the (/)^ 
theory, at the moment of time reflection symmetry, t = 0: 



S 



=0) 



=0) 



--0) 



3 
1 
9 

1 
9 

1 
9 

2 

27 



195 



35 



65 



SZ+ 10 



Y 



2275 
64 



dS 
dp 



S^Z 



26D{D 



325 dS dZ 
4 dp dp 

1) fdS_ 
dp 



8 „9 

+ 3^ - 


125 ^4 
T2"* 






1503^ 


-fA 


'dsV 


16 


ydp) 


4225^2 


8125 


S^Z + 


128 


48 



^ - 245* - 



595 
"96" 



53 



11285 

384 



(48) 
(49) 
(50) 

(51) 

(52) 



6589 



SZ- 



9223 



52{D - 1] 



^dS 
dp 



92 
~9' 



32 

35417 

" 288 



dS 
dp 



dSV 



dp J 



21467, 
144 



(53) 



These expressions are presented because they are needed to provide good initial data for numerical time evolution 
simulations. They will be actually used in the next Section llVI Since the range of 0^,^"'^^ increases very much with 

k, we depict the product s''4>^^~^^ for some chosen values of e on Figs. [51 El [7] and [51 On these four figures Sg and 
Si are depicted in spatial dimensions D = 2 and D = 3. In a given order in the expansion, the value of e which 
brings the contribution of lower order terms approximately to the same order of magnitude will be used as an upper 
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estimate for the range of e, below which our expansion can still be expected to yield an acceptable approximation. 
Generally speaking, in an asymptotic expansion for a given (small) value of the expansion parameter one can only 
sum terms up to such an order until which all terms decrease. A simple comparison of a QB corresponding to the 

D=3 , Si 

0.5 
0.4 
0.3 
0.2 

o~ 

II 

^ 0.1 


-0.1 
-0.2 
-0.3 

0.5 1 1.5 2 2.5 

P 

FIG. 6: e*^^^^'''' corresponding to the solution 5*1, with one node in the three dimensional case. Because of the sharp increase 
of the higher order functions the value of e was chosen to be a much smaller as for Sq. 

fundamental solution, Sq, with another one corresponding to a solution with a single node, 5*1, (compare Figures [5] 
and[6l), makes one to guess that oscillons containing fundamental QB's are likely to have better stability properties 
and longer lifetimes, than those containing QB's based on solutions with nodes, at least in Z? = 3. Our numerical 
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simulations show that this is indeed the case (see Section |IV|) . In two dimensions the difference between the longevity 
and stability properties of oscillons containing QB's corresponding to solutions of the S equation (PT|) with nodes is 
much less pronounced than in three spatial dimensions. It is also apparent that the e expansion is valid for significantly 
larger values of e in the two dimensional case than in the three dimensional one. 





In order to get an independent check on the validity of the small amplitude expansion, we have compared some 
QB's up to order 6, with time periodic QB's obtained previously by solving the NLWE ^ directly by Fourier mode 
decomposition flo|. On Figures [51 and [TUl we depict (j) computed to various orders in the e-expansion, and also the QB 
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obtained in Ref. [10| by Fourier mode decomposition using very precise spectral methods provided by the LORENE 
Hbrary [3l|. The chosen frequencies correspond to two states investigated in detail in |10|]. Note that in Ref. [lOj 
the interaction potential had a different scale, and the resulting threshold frequency was \f2 as opposed to the value 
1 in the present paper. The periodic quasi-breather solutions chosen from Ref. [l3| have frequencies a) = 1.412033 
and w = 1.398665 (see Figs 5. and 19. in [l3|)- In the present conventions these values correspond to frequencies 
w = 0.9984581 and uj = 0.9890055. The corresponding values of e are: e = 0.05551039 and e = 0.1478787. It can 
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FIG. 9: Comparison of the value of the field at r = obtained by the e expansion to a very precise value obtained by a high 
order Fourier mode decomposition for uj = 0.9890055, i.e. for e = 0.1478787 in D = 3. For such a large value of e the first 
order approximation gives a remarkably good estimate in the neighbourhood of the central region. The contributions of the 
higher order approximations make this agreement increasingly worse in the neighbourhood of the origin. Farther away from 
the center, however, the third order approximation gives the best, although not a very precise, result. 



be seen on Figures 111 [TU] that for higher values of e, only the leading term or eventually the first two orders in the e 
expansion give meaningful results. The relatively big error in these approximation cannot be decreased because of the 
asymptotic nature of the expansion. As e gets smaller and the frequency w gets closer to the basis frequency 1, more 
and more higher order terms in the e expansion can be used, and then the error also decreases significantly. It can be 
seen on Fig. [TO] that for the smaller value e although the second, third and fifth order expansion gives an improvement 
on the lower order values, the fourth and sixth order expansion turns out to be less precise than the third and fifth 
order expressions. We think that this is related to the fact that the signature of the various order contributions 
changes in pairs, i.e. 0^^"°^ is positive at the center r = for fc = 1,2,5,6... and negative for k = 3,4,7,8..., as can 
be seen on Fig. [51 This alternating improving and not improving behaviour for odd and even orders happens only for 
intermediate values of e. For even smaller e, e.g. for e = 0.01, the error decreases monotonically when increasing the 
order of the expansion. 

Let us now come back to the e-expansion of the energy computed in subsection III Dl From our results it is 
not difficult to calculate numerically Eq and Ei in D = 2 and D = 3. We find that the first two terms of the time 
averaged energy in the e-expansion are given as: 

3.9003 + 26. 9618e^ for Z) = 2 , « 6.29908/e + 264.262e , for D = 3 . (54) 

This simple estimate Eq. (j54p gives for the minimal value of e^n ~ 0.15428, which is unfortunately already too large to 
be trusted. Nevertheless it can still be accepted as the indication that such a minimal value, lu^ exists. For spherically 
symmetric oscillons in D = 3 it has been found that ~ 0.9659 [l3|- This value of Wm corresponds to ~ 0.2588 
which value is way too large for us. 
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FIG. 10: Comparison of various order e expansions of at r = to the value obtained by Fourier mode decomposition for 
w = 0.9984581 i.e. for e = 0.05551039 in the D = 3 case. 



IV. TIME EVOLUTION 



The precision and applicability of the e expansion can be checked by using the field value obtained by the expansion 
as initial data for a numerical time evolution code applied for our spherically symmetric scalar field system. The field 
value given by Eqs. - ([55]) and ([5]) at r = t = has been used as initial data for various e values in I? = 2 and 
D — S spatial dimensions. 

The applied numerical evolution code is a slightly modified version of the fourth order method of line code used in [l3| 
for studying oscillons and developed in [35| for the study of spherically symmetric magnetic monopole configurations. 
The spatial grid is chosen to be uniform in the compactified radial coordinate R defined by 

2R 

where k is a constant which may be chosen differently though for each choice of initial data. The whole range 
< r < oo of the physical radial coordinate r is mapped to the interval < i? < 1, avoiding the need for explicitly 
describing boundary conditions at some large but finite radius. Since the characteristic size of the obtained oscillon 
states is inversely proportional to e we chose k to be proportional to e, keeping the oscillon occupying approximately 
the same region of the R coordinate range. For the actual calculations we used k = 5e. 
Using the radial coordinate R the field equation ^ takes the form 



2\3 



n'jl-R') 

2(l + i?2) 



2(l + i?2) 



Introducing the new variables 



ft 



i?(3 + i?2 

(l+i?2)2 



2R 



-£/'(</)), 



\R 



(56) 



(57) 
(58) 



16 



the problem can be interpreted as a system of first order differential equations comprising ([57]) . (|56p. and = (j)t^R 
for the three variables 0, (pt and (jjn. If equation (|58p holds at i = it is preserved by the evolution equations, 
thereby it can be considered as a constraint. The third term inside the bracket on the right hand side of (|56p cannot 
be directly evaluated numerically at the center R = 0. At the grid point corresponding to the center this term is 
calculated using the identity 

r^O Jx r— ►O 

Since the evolution of the initial data provided by the e expansion show different characteristics in two and three 
dimensions we discuss these cases in different subsections. 



A. D = 3 



In the three dimensional case, there are two different types of oscillons, a stable and an unstable type. For the 
frequency range w < ~ 0.967, (i.e. for e > 0.255), oscillons are essentially stable. They slowly radiate energy 
while their frequency, uj{t), increases towards a critical frequency lOq- When they reach the critical frequency these 
oscillons quickly disintegrate. Oscillons with ut > Uc have one unstable decay mode, which can be suppressed by fine 
tuning the initial data. Close to the critical value of the parameter in the initial data there can be two types of decay 
mechanisms. One with a uniform outwards motion of the energy, and another through a temporary collapse to a 
small central region (see Fig. 4 of llOl ) . In order to find such oscillons a very precise fine tuning of the initial data 
is necessary. For example in Ref. |10j this fine tuning corresponded to the classical bisection procedure between two 
values of a suitable parameter in the initial data yielding the two different decay modes. This way one obtains very 
long living oscillon states. The frequency of these unstable oscillons decreases slowly towards uj{t) — > uic- 

Remarkably the energy of QB's as a function of their frequency exhibits a minimum at a; « Wc, (see Fig. 3 in 0], 
Fig. 17 of [13] and Fig. 4 of [2§|)- Therefore it is natural to assume that the two types of oscillons are also distinguished 
by the same behaviour of their energy as function of their oscillon frequency. 

Unfortunately, for values of e > 0.25 initial data obtained by the e expansion are well outside the domain of validity 
of the expansion. Using such initial data gives decaying states which are unrelated to the stable oscillons with the 
intended frequency. On the other hand the e expansion yields good initial data for small amplitude unstable oscillons. 
As a matter of fact the initial data obtained this way makes a fine-tuning procedure of the initial data unnecessary. 
For sufficiently small values of e < 0.1 the first few terms (at least to order 0{e^)) of the series expansion - ([55]) 
already yield sufficiently good initial data which evolve directly into long living oscillon states. 

For the three dimensional case we present the results of the time evolution of the initial data obtained up to order 
six by the e expansion method for two choices of e. First we consider initial data obtained from the basic solution Sq 
of pT|) without nodes. As we will see shortly, those with nodes provide initial data that leads to states of significantly 
shorter lifetimes. Figure [TT] shows the upper envelope of the central value of cj) for the initial data belonging to 
e = 0.05551039. The decay method changes with the order of the initial data. The amplitude peak on the evolution 
of the third and fourth order initial data reflects the collapsing decay mode of the oscillon state. Although, in general, 
lifetimes get longer for higher order approximations, this increase is not monotonic. In accordance with Fig. llOi 
approximations of order 4 and 6 do not bring any improvement on the the functions of order 3 and 5. The time 
evolution of the frequency of the oscillations is plotted on Fig. [T^ It can be seen that for this relatively high value of 
e the first and second order approximation yields a shorter living state quite different from the expected oscillon state 
with frequency ui — 0.9984581. On Figure [13] the evolution of a higher frequency initial data with e — 0.01 is shown. 
It can be seen that initial data with smaller e provide evolutions with significantly longer lifetimes. Fig. [TJ] shows 
the initial stage of the evolution. It can be seen that although for a short time the error of the solution decreases 
monotonically with the order of the e expansion, for longer time intervals the fourth and sixth order approximations 
do not improve on the previous order expansions. 
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FIG. 11: Time evolution of the initial data obtained by the e expansion method up to order six for e — 0.05551039 in three 
spatial dimensions. 
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FIG. 12: Oscillation frequency as a function of time for the states shown on Fig. 1111 




FIG. 13: Time evolution of the frequency of the oscillations of 4> evolving from initial data with e = 0.01 in the D = 3 case. 
Since the frequency is very close to one, the value a; — 1 is plotted instead of ui. 
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FIG. 14: Beginning of the time evolution simulations shown on Fig. 1131 For this shorter time interval the difference from the 
expected solution with u) = 0.99994999875 decreases monotonically with the order of the e approximation. 
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Figure [15] shows the evolution of various order initial data with the same e as on the previous two figures but 
obtained by using the solution 5*1 of (PTjl with one node. These are localized, although big size, high energy states. 




FIG. 15: Time evolution of the initial data obtained by the e expansion up to order four using the solution 5i of (|2ip . Space 
is three dimensional and e — 0.01. The top graph shows the upper envelope of (j) at r — 0, while the lower graph plots the 
evolution of the oscillation frequency. It can be seen that the first and second order initial data still evolve to states very 
different from the one with the expected frequency. 



The time dependence of these states is rather complex; there is an interior part (r < 50) which after an initial time 
interval (w 1200) shows a complicated time dependence, with large amplitude variations. The time dependence of 
these states is close to being periodic for recurrent time intervals. In the initial stages these initial data still evolve 
close to a periodic configuration with the expected frequency, although it can stay near this state much shorter time 
than the evolution obtained using Sq. Decreasing the value of e increases the lifetime of these states as well, but they 
still remain less stable and shorter living than the basic states obtained by using Sq. 

On Figure [16] the value of as a function of the radial coordinate r for subsequent time slices is plotted during a 
half period of oscillation. For comparison, the corresponding configuration obtained from initial data generated with 
5*0 is also presented. In case of Si initial data the value of (fi remains very close to zero at r = 51.3. The energy 
density £ also remains very small at this radius, as can be seen on Fig.[T71 The total energy of the Sq configuration 
is E = 632.536, while the Si configuration contains significantly more energy, E ~ 4061.88. 
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FIG. 16: Time evolution of (j) during a half period of oscillation, plotted as a function of the radial distance r. The upper plot 
corresponds to initial data obtained from So, the lower from Si. The absolute value of at r = 51.3 remains below 10~^. 
Both plots contain lines corresponding to uniform time steps between the first maximum and minimum after t — 300, although 
the plot would remain very similar in a large time interval. The initial data in the So case was generated by a sixth order e 
expansion, while a third order expansion was used in the Si case. 
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FIG. 17: Energy density £ of oscillon configurations evolved from So and Si initial data, using e = 0.01. The energy density 
is plotted on the same uniformly placed moments of time as on Fig. 1161 The minimum of £ in the Si case remains below 10~^ 
near r = 51.3. 
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B. D = 2 



In two spatial dimensions there seems to be a single type of oscillon, which is stable. Once formed, oscillons in 
D = 2 are not observed to disintegrate. Their energy is a monotonically decreasing function of the frequency. This 
explains the observation that all oscillon states are stable. When they slowly emit energy by radiation they gradually 
evolve through oscillon states with increasing frequency, U!{t) — > 1. Oscillons evolve from a wide range of initial data, 
by shedding most of the surplus energy quickly during an initial state. However, in general, a slow periodic change 
can be seen on the amplitude and on the frequency of the oscillations, indicating a breathing type oscillation of the 
oscillon as a whole. The amplitude of this low frequency ringing depends on how closely the initial data approaches 
a given pure oscillon state. 

Since, similarly to the D — 3 case, evolutions from initial data obtained using solutions S of ((2T|) with nodes 
produce less stable and shorter living states, in the following we present only numerical simulations corresponding 
to the nodeless solution Sq. In two dimensions the validity domain of the expansion extends to significantly higher 
values of e. On Fig. [18] we present time evolution results for e = 0.218632 corresponding to u — 0.9758073. It can 
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FIG. 18: Time evolution of the initial data produced by various order e expansions for e = 0.218632 in the two dimensional 
case. The top graph shows the upper envelope of at r = 0, while the bottom graph the time dependence of the frequency. 
The evolution of the third and fourth order initial data clearly corresponds to some oscillon state deformed by a low frequency 
ringing mode. The evolution of the fifth order initial data corresponds to a similar state with extreme high amplitude ringing. 
The first, second and sixth order initial data evolve into decaying modes, indicating a large error at these orders of the e 
expansion. 

be seen that at this high e value the approximation improves up to order tree in the expansion, and then starts to 
deteriorate in line with the asymptotic nature of the expansion. 

On Fig. [in] the e = 0.05551039 case is presented. This value of e with the corresponding frequency u> — 0.9984581 
has been also studied in D = 3 dimensions. In this intermediate frequency case the first and second order initial 
data still lead to decaying evolutions. Higher order expansions tend to give improving approximations of an oscillon 
state ringing with a low frequency. Increasing the order of the expansion gives smaller amplitude ringings, thereby 
approaching a pure, very closely periodic, oscillon state. Similarly to the three dimensional case, initial data of order 
four and six yield larger error than order three and five, which, however, does not change the overall improving 
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FIG. 19: The frequency as a function of time for the evolution of initial data with e = 0.05551039 in case of D = 2. 
tendency of the initial data. 

The average amphtude and average frequency of the presented states show extremely little change even for much 
longer time periods than the ones presented on the figures. The amplitude of the low frequency ringing decreases very 
slowly too. To determine the rate of change of the amplitude and the corresponding slow energy loss by radiation 
would require very high resolution numerical runs requiring excessive processor time. 

V. CONCLUSIONS 

Small amplitude oscillons represent an important subset of time dependent long- living lumps. We have shown that 
they can be very well approximated by an asymptotic series of localized, time-periodic breather-like objects (quasi- 
breathers). We have developed a general framework to derive the asymptotic series expansion of small amplitude 
quasi-breathers, in D spatial dimensions in general scalar theories. We have derived a 2nd order elliptic PDE with a 
cubic non-linearity, universal for scalar models, which determines these quasi-breathers. Our numerical investigations 
in ^*-theories show that the small amplitude quasi-breathers obtained by the asymptotic expansion, provide excellent 
initial data for long-living oscillons in Z) = 2 and D = 3. We have found that small amplitude QB's do not exist for 
D> 4. 
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